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ABSTRACT 

We use a combined 120ks Chandra exposure to analyze X-ray edges produced by non-hydrostatic gas mo- 
tions (sloshing) from galaxy collisions, and cavities formed by AGN activity. Evidence for gas sloshing is seen 
in the spiral morphology and multiple cold front edges in NGC 5846's X-ray surface brightness distribution, 
while lack of spiral structure in the temperature map suggests the perturbing interaction was not in the plane 
of the sky. Density and spectral modeling across the edges indicate that the relative motion of gas in the cold 
fronts is at most transonic. Evidence for AGN activity is seen in two inner bubbles at 0.6 kpc, filled with 
5 GHz and L5GHz radio plasma and coincident with Ha emission, and in a ghost bubble at 5.2 kpc west of 
NGC 5846's nucleus. The outburst energy and ages for the inner (ghost) bubbles are ^ 10^^ ergs and ^ 2 Myr 
('^ 5 X 10^^ ergs and 12Myr),respectively, implying an AGN duty cycle of 10 Myr. The inner bubble rims are 
threaded with 9 knots, whose total 0.5-2keV X-ray luminosity is 0.3 x lO'^^ergs"^ , a factor ^ 2-3 less than 
that of the surrounding rims, and 0.7 keV mean temperature is indistinguishable from that of the rims. We sug- 
gest that the knots may be transient clouds heated by the recent passage of a shock from the last AGN outburst. 
We also observe gas stripping from a cE galaxy, NGC5846A, in a 0.5 kpc long (^ IQt'Mq) hot gas tail, as it 
falls towards NGC 5846. 

Subject headings: galaxies: clusters: general - galaxies:individual (NGC 5846) - galaxies: intergalactic 
medium - X-rays: galaxies 



1. INTRODUCTION 

Over the past decade, deep X-ray observations of signifi- 
cant substructure in the hot gas near the centers of rich clus- 
ters have revolutionized our understanding of the evolution of 
clusters and the galaxies at their cores. Sharp surface bright- 
ness discontinuities (edges) are found where the increase in 
density across the discontinuity is accompanied by a de- 
crease in temperature, opposite to that expected from shocks. 
These 'cold fronts' provide evidence for an ongoing subclus- 
ter merger or non-hydrostatic motions of the gas ('sloshing') 
from a prior disturbance in the core (see, e.g. the review by 
Markevitch & Vikhlinin 2007). Furthermore, episodic out- 
bursts from the active galactic nucleus (AGN) at the centers 
of the dominant galaxies in groups and clusters, are revealed 
by outflow shocks and X-ray cavities (deficits in X-ray sur- 
face brightness) in the intra-cluster medium (ICM) (see, e.g. 
Fabian et al. 2003 for the core of the Perseus cluster, Forman 
et al. 2005 for M87 in the Virgo Cluster). These X-ray cavi- 
ties correspond to regions where the hot cluster gas has been 
evacuated by jets from AGN outbursts and the interior of the 
cavity (bubble) is filled with radio plasma. After an initial mo- 
mentum driven formation stage, the 'bubble' of radio plasma 
rises buoyantly in the group/cluster atmosphere. This feed- 
back between the central AGN and its environment may be 
key to resolving the cooling flow problem in galaxy clusters, 
and and may be an essential ingredient in the emerging mod- 
els for the coevolution of galaxies and their central supermas- 
sive black holes (e.g. Best et al. 2006; Croton et al. 2006). 
While these phenomena have been studied in high angular 
resolution, deep X-ray observations for many rich clusters 
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(e.g. McNamara et al. 2000; Nulsen et al. 2002; Birzan 
et al. 2004; Nulsen et al. 2005; McNamara et al. 2005; Allen 
et al. 2006), their impact on structure formation in the much 
more common cores of galaxy groups has received far less 
attention (however, see e.g. Machacek et al. 2010 for the 
Telescopium group; Gitti et al. 2010 for HCG62; Randall 
et al. 201 1 for the NGC 5813 group; David et al. 201 1 for the 
NGC5044 group). Yet, it is in groups at high redshift where 
galaxy evolution is most rapid (Cooper et al. 2006), gas-rich 
mergers are more common, and, since the gravitational po- 
tential is shallower, outbursts from supermassive black holes 
(SMBHs) may have a greater effect on the surrounding gas 
than in their richer cluster cousins. 

In this paper we present results from the combined ~ 
120 ks Chandra observation of the NGC 5846 group. The 
Digital Sky Survey (DSS) image is shown in Figure [T] 
Since the NGC 5846 group is only 24.9^2 2 Mpc distant 
(Tonry et al. 2001), this exposure is one of the deepest 
high angular resolution X-ray observations of any nearby 
group of galaxies. The NGC 5846 group consists of ^ 20 
bright (Mb < -14.9) galaxy members with a group veloc- 
ity dispersion of 386!|^kms-' (Zabludoff & Mulchaey 1998; 
Mulchaey et al. 2003). Previous X-ray studies using the 
Einstein X-ray Observatory (Forman, Jones & Tucker 1985; 
Biermann, ICronberg & Schmutzler 1989) and using ROSAT 
and ASCA (Finoguenov et al. 1999 and Trinchieri et al. 
1997), show ^ 1 keV diffuse group X-ray emission extend- 
ing to ^ 200 kpc (^ 28 f7) from the central, brightest group 
galaxy NGC 5846 and a total X-ray luminosity within 200 kpc 
of ^ lO^^^ergs'^Mulchaey et al. 2003). From combined 
optical and X-ray studies using ROSAT, ASCA and XMM- 
Newton data, Mahdavi, Trentham & Tulley (2005) found that 
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Fig. 1 . — DSS image Bj band image of the NGC 5846 group taken with the 
UK Schmidt Telescope on 1980 March 20. North is up and East is to the left. 



the NGC 5846 group may be part of a larger (8 x lO'^M©) 
structure, sharing a common 1.6Mpc diameter dark matter 
halo with the NGC 5813 group located 600 kpc away in pro- 
jection, and that group membership is overwhelmingly dom- 
inated by dwarf elliptical galaxies (80%), similar to that seen 
in the Virgo cluster There is, however, no evidence for inter- 
action between the NGC 5846 and NGC 5813 groups (Randall 
etal. 2011). 

NGC 5846 is a giant EO elliptical galaxy with an effec- 
tive radius of 11.6 kpc (Bender et al. 1992) and an abso- 
lute V magnitude of -22.6 (Faber et al. 1989), compara- 
ble to M87 (My = -22.5), the giant elliptical galaxy in the 
Virgo Cluster However, in contrast to M87 whose central 
velocity dispersion suggests a black hole mass of at least 
~ 3 X 10^ Mq, NGC5846's central velocity dispersion of 
278 kms"' corresponds to a supermassive black hole mass of 
3.5 X IO'^Mq, an order of magnitude smaller (Tremaine et al. 
2002; Gebhardt & Thomas 2009). Radio emission from the 
nuclear region of NGC 5846 is complex, with a central com- 
pact, possibly variable, steep spectrum (a = 0.8) radio source 
with radio fluxes of 1.5 mJy (2.8 mJy) at 5 GHz (2.3 GHz), 
likely associated with an AGN. Three additional radio com- 
ponents are seen in 2.3, 5, and 15 GHz VLBA maps roughly 
aligned in the north-south direction with the nuclear emis- 
sion (Filho et al. 2004). These regions, extended on milli- 
arcsecond scales, may either be compact supernova remnants 
or emission from a jet (Filho et al. 2004). XMM-Newton 
observations suggest at least 3 surface brightness discontinu- 
ities (edges) in the gas surrounding NGC 5846 (Finoguenov 
et al. 2006). Such discontinuities indicate gas motions. These 
may be residual gas 'sloshing' from a previous galaxy merger 
or from a high velocity galaxy encounter with the large, HI 
rich, barred spiral galaxy NGC 5850, located 10.'3 (7 1.5 kpc 
in projection) to the east (see Fig. [T]). The large shift of HI 
gas to the west and northwest in NGC 5850, in the direction of 
NGC 5846, and NGC 5850's distorted stellar ai-ms and spurs, 
also to the west, suggest a recent ( < 200 Myr) tidal interaction 
between the two galaxies (Higdon et al. 1998). 

Non-hydrostatic gas motions also could be induced by 
AGN activity. A previous short (30 ks) Chandra X-ray obser- 
vation of NGC 5846 shows bright, knotty -rimmed X-ray cavi- 
ties associated with radio lobes in NGC 5846's central ~ 2 kpc 
(Trinchieri & Goudfrooij 2002; Allen et al. 2006). Extensive 
filamentary Ha + [Nil] line emission and dust also are found 
in the central region of the galaxy, highly correlated with 
the X-ray emission (Goudfrooij & Trinchieri 1998). Spitzer 
MIPS observations of NGC 5846 show extended 70/im dust 
emission to a FWO.IM ^ 5 kpc, suggesting that dust, origi- 



nating from stellar mass loss from evolved AGE stars, accu- 
mulated in the cental ^ 1 kpc of that galaxy and was subse- 
quently transported buoyantly to these larger radii by AGN 
activity (Temi et al. 2007). 

In this paper we use the combined 120 ks Chandra exposure 
to study in detail the X-ray cavities and surface brightness 
edges that chronicle AGN activity in NGC 5846 and the rela- 
tive motions between the galaxy ISM and group gas. In ^we 
summarize our observational data reduction and analysis tech- 
niques. In |3] we analyze the large scale surface brightness 
and temperature structure in the core of the NGC 5846 group, 
discuss evidence for gas sloshing, and model the density, 
temperature and pressures across observed surface brightness 
edges. In ^ we focus on the properties of the system of 
AGN outburst cavities near the center of NGC 5846. In ^ 
we also present the first observation of gas stripping from a 
compact elliptical galaxy, NGC 5846A, during its' infall to- 
wards NGC 5846. Our conclusions are summarized in ^ 
Unless otherwise indicated, quoted uncertainties are 90% CL 
for spectra and 1 a for photon counts. For the standard A dom- 
inated cold dark matter cosmology (Hq = 71 kms"' Mpc"' , 
Jl,,, = 0.27, JIa = 0.73) and assuming NGC 5846 is at rest at 
the center of the group potential, the luminosity distance to 
NGC 5846 is 24.2 Mpic, and 1" = 0.116kpc (Wright 2006). 
This is consistent within uncertainties with the luminosity dis- 
tance determined from surface brightness fluctuations (Tonry 
2001). 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 5846 was observed with the Chandra X-ray Observa- 
tory for 30313 s on 2000 January 24 (obsid 788) in FAINT 
mode with ACIS-S at aimpoint, and for 90367 s on 2007 June 
12 (obsid 7923) in VFAINT mode with ACIS-I at aimpoint 
for a total of ^ 120 ks. The data were initially filtered to re- 
ject bad pixel grades (1, 5, 7) and data that fell on hot pixels. 
Data flagged in obsid 7923 (VFAINT mode) as having ex- 
cessive counts in the border pixels surrounding event islands 
were also removed to optimize signal to noise at energies be- 
low 1 keV. We then reprocessed these data with CIA04.2 us- 
ing standard tools and calibrations provided by the Chanda 
X-ray Center, applying corrections for the charge transfer in- 
efficiency on the ACIS CCDs, the time-dependent build-up of 
contaminant on the optical filter, and the secular drift (tgain') 
of the average pulse-height amplitude for photons of fixed en- 
ergy. Periods of anomalously high particle count rates (flares) 
were excised using lc_clean and a 3cr clipping algorithm, re- 
sulting in useful exposures of 23205 s and 89509 s for obsid 
778 and obsid 7923, respectively. 

Backgrounds for the imaging and spectral analyses, where 
indicated, were drawn from the source free background data 
set^, i.e. period C ACIS-S ( period D ACIS-I) for OBSID 788 
(7923), appropriate for our observation and instrument con- 
figuration, using the CI AO acis_bkgrnd_lookup tool. These 
data were reprocessed using the same gain files as used for 
the source data, reprojected to the aimpoint of the observa- 
tion, and then normalized by comparing count rates in the 
9.0- 12.0keV energy band, where particle background domi- 
nates. This resulted in additional normalization factors of 0.94 
and 1.44 for obsid 788 and obsid 7923, respectively. Since 
NGC 5846 lies close to the Northern Spur, there is an addi- 
tional soft Galactic background component not accounted for 

' see, http://cxc.harvard.edu/contrib/alexey/tgain/tgain .htmll 
^ see http://cxc.harvard.edu/contrib/maxim/acisbg 
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Fig. 2. — Chandra background-subtracted, exposure-corrected, coadded 
0.75-2.0 keV image of diffuse gas in the outer regions of NGC 5846 showing 
the main X-ray features outside 2kpc of the nucleus. Ipixel = l"0 X l"0 and 
the image has been smoothed with a 2"0 gaussian kernel to highlight faint 
features. North is up and east is to the left. 



in the blank sky background data. This Galactic component 
is expected to be uniform over the combined ACIS fields- 
of-view (FOVs). We model this background using the flux 
measured in a 63 "55 radius circular region on CCD6 cen- 
tered at (RA, Dec = 15''06'"56.391'', -l-01°20'33"660, a dis- 
tance of 17fl6 from NGC 5846 where group emission is ex- 
pected to be small. This soft Galactic background component 
was also subtracted from the images analysed in ^ Instru- 
ment and exposure maps were generated in 0.25 keV pass- 
band intervals from 0.5 to 4.0 keV and 2keV intervals from 
4.0- lO.OkeV for both data sets using standard CIAO tools. 
Background-subtracted, exposure corrected images were con- 
structed in various bandpasses and the data sets coadded for 
further analysis. 

The CIAO tool wavdetect was run on the coadded, 
background-subtracted 0.5 - 4.0 keV image to identify point 
sources, and the source regions were adjusted to eliminate 
overlapping regions. An ellipse slightly larger than the source 
region was generated as a background region for each source. 
Sources outside the central 2.3 kpc (20") of NGC 5846, i.e. 
those not associated with knots on the inner bubble rims 
or NGC5846's nucleus, were excised from the data. The 
point source regions in the exposure-corrected images were 
then filled with events sampled from the corresponding back- 
ground region around each source to produce the images of 
diffuse gas shown in ^ 

Spectra were extracted from the point source cleaned event 
files using standard CIAO tools and fit with XSpec 1 1 .3.0. Lo- 
cal background regions drawn from the ACIS-I detector are 
used, when possible, to mitigate the effects of the soft Galac- 
tic background. The spectra from both observations were fit 
jointly, when such a fit reduced uncertainties in the model pa- 
rameters. Otherwise, we present only the long 89 ks exposure 
(OBSID 7923) for our spectral analysis, because of its higher 
number of source counts. 

3. LARGE SCALE STRUCTURES 

In Figure |2] we show the 0.75- 2 keV X-ray background 
subtracted, exposure corrected, coadded image of diffuse gas 
in and surrounding NGC 5846. The image has been binned 



by 2 pixels into ^1" bins, and smoothed with a 2."0 Gaus- 
sian kernel to highlight faint emission features outside the 
nuclear region of NGC 5846. The surface brightness distri- 
bution is brighter (gas denser) from ^ 4 -20 kpc northeast of 
the nucleus than at similar distances to the southwest. Surface 
brightness discontinuities (edges) are found ^ 2.'7 19 kpc) 
to the northeast and southwest of NGC5846's nucleus, with 
possible inner edges to the southwest at '-^ lf5 ( ^ 10.5 kpc) 
and to the northeast at ^ 1' ('^ 6.7 kpc). In 33.3. ll we find that 
the ratio of gas density across these edges exceeds 1 at > 90% 
confidence level in all cases, confirming that these disconti- 
nuities are significant. We also observe a tail-like feature of 
enhanced emission (NW extension) at 13(t significance in the 
flux over expected background, that extends 4.'1 (28.5 kpc) to 
the northwest suggesting a spiral pattern. A roughly spherical 
deficit in X-ray surface brightness is found, at an 8(t signifi- 
cance in flux compared to neighboring regions, located ^ 45" 
(5.2 kpc) west of the nucleus. In contrast to the inner X-ray 
cavities observed in previous studies (Trinchieri & Goudfrooij 
2002; Allen et al. 2006), this X-ray cavity is not found coin- 
cident with radio synchrotron emission (see ^and Fig. [TTI ). 
Thus it is a 'ghost bubble' where the electron population in 
the radio plasma, expected to fill the bubble, has aged such 
that the synchrotron emission has moved to lower frequencies 
that are more difficult to observe (Birzan et al. 2004). Diffuse 
X-ray emission is also observed from the neighboring galaxy 
NGC 5845. 

3.1. Gas Densities: Evidence for Sloshing 

To investigate these asymmetries in more detail, we use 
SHERPA to fit a single component /3-model to the mean 
emission from NGC 5846 and then study deviations from this 
mean. In the upper panel of Figure[3j we show the 0.5-2 keV 
azimuthally, averaged radial surface brightness profile, cen- 
tered on the nucleus of NGC 5846 superposed with the best-fit 
single /3-model (solid line). We find that within 4 kpc the sur- 
face brightness profile is best fit by a single /3-model with pa- 
rameters (3 = 0.452, core radius rc = 1.08 kpc (9"3), and cen- 
tral surface brightness io = 8.95 x 10"^ photons s"' cm"^ , cor- 
responding to a central electron density of no = 0.091 cm"^ . 
While j3 is in reasonable agreement with previous ROSAT re- 
sults ifi ^ 0.5), we find a factor ^ 1.6-2 smaller core radius 
and infer a factor 2 lower mean central electron density com- 
pared to previous results (Trinchieri et al. 1997; Osmond & 
Ponman 2004). These discrepancies are likely a reflection of 
the different angular resolution of the detectors and the fail- 
ure of a single /3-model to describe the asymmetric surface 
brightness distribution. 

In Figure |4] we present the fractional difference map be- 
tween the observed data and the best-fit 1 -dimensional /3- 
model for the group gas emission (i.e. left panel of Fig. 
[3]). We see clearly the characteristic spiral pattern in surface 
brightness (and thus density) generated by 'sloshing'. Gas, 
initially displaced by the gravitational attraction and/or hydro- 
dynamical forces from the off-axis encounter, is set into oscil- 
lation as the ram-pressure suddenly decreases and the gravi- 
tational potential energy and angular momentum imparted to 
the galaxy gas by the close encounter causes the higher den- 
sity core gas to rebound and overshoot the center of the group 
dark matter potential (black cross in Fig. |4] east of the X- 
ray 'ghost' cavity) and also spiral outward (see, e.g. simula- 
tions by Ascasibar & Markevitch 2006; ZuHone et al. 2010). 
The faint northwest extension, identified in Figure |2] and seen 
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Fig. 3. — (upper) Mean radial surface brightne.ss profile of NGC5846 in 
the 0.5 — 2.0keV energy band. The black curve is the best fit single /3- 
model with /3 = 0.452, core radius n- = 1.08 kpc {9'.'3) and central surface 
brightness .so = 8.95 X 10"^ photons s"' cm"^ . 1 pixel = 1"0 X I'.'O. (lower) 
0.5 — 2.0keV surface brightness radial profiles measured from 123° to 180° 
to the northeast (black squares) and from 268° to 338° to the southwest (gray 
circles) compared to the mean beta model shown in the left panel (dashed 
line). All angles are measured counterclockwise from west. 



clearly in Figure|4] is the tail of this spiral structure, as denser 
core gas is pushed to larger radii in response to the encounter. 
This tail-like structure is similar to that first found in the X- 
ray surface brightness map of the NGC 5098 group (Randall 
et al. 2009). In the lower panel of Figure |4] we show the 
simulated X-ray emission from a cluster merger adapted from 
Figure 19 of Ascasibar & Markevitch (2006), showing sim- 
ilar spiral density structure and a tail-like density enhance- 
ment caused by the non-hydrostatic bulk motions (sloshing) 
of gas in the cluster core. Since the simulation was for a 
more massive cluster merger we do not expect the scales to 
be the same. However, the similarity of our observation of 
sloshing in groups with features in simulations of galaxy clus- 
ter mergers demonstrates that such surface brightness features 
are common to off-axis galaxy, group and cluster encounters 
over a wide range of mass scales (see, e.g.. Fig. 1 and Fig. 9 
in Johnson et al. 2010 for galaxy cluster Abell 1644; Fabian 
et al. 2011 for the Perseus cluster; Randall et al. 2009 for 
the NGC 5098 galaxy group; Machacek et al. 2010 for the 
NGC 6868 group). 

However, as also seen in Figure 19 of Acasibar & Marke- 
vitch (2006), the projected surface brightness and correspond- 
ing temperature features are highly dependent on the viewing 
angle between the plane of the encounter and the observer. If 
one ignores the NGC 5846 ghost bubble and the faint north- 
western surface brightness tail, the image of NGC 5846 in Fig. 
|2] shows two ellipsoidal edges to the southwest and, although 




Fig. 4. — (upper) Fractional difference map ((Im — Mod) /Mod) between 
the 0.5 — 2keV X-ray surface brightness image of the NGC 5846 group (Im) 
and the best-fit 1-dimensional /3-model (Mod) for the data (see Fig. [3). 
Both Im and Mod were smoothed with a T.'5 Gaussian kernel. The nu- 
cleus of NGC 5846 is denoted by a black cross and the ^ 30" (3.4 kpc) di- 
ameter X-ray 'ghost' cavity is outlined by the circular contour to its west. 
1 arcmin = 7 kpc and the panel size is 90.5 kpc. Note that the linear sequence 
of blobs in the spiral tail are artifacts due to the ACIS-I chip gap. (lower) The 
predicted X-ray emission from a simulated galaxy cluster merger, that shows 
similar structure caused by 'sloshing'. The image was adapted from Fig. 19, 
Ascasibar & Markevitch (2006), and the panel size is 250 kpc. 



the inner features are less pronounced, at least two edges to 
the northeast. Such multiple surface brightness edges are also 
characteristic of gas sloshing (Markevitch et al. 2001) and 
are similar to the y projection of simulated cluster mergers 
in the right panel of Figure 19 from Ascasibar & Markevitch 
(2006). The orientation of the plane of the encounter, that 
caused gas sloshing in NGC 5846, with respect to the ob- 
server's line of sight likely lies between the z and y simulation 
projections, such that the sloshing spiral is tilted with respect 
to the plane of the sky. This is consistent with the identifica- 
tion of the perturber as the disturbed spiral galaxy NGC 5850, 
whose relative radial velocity with respect to NGC 5846 is 
large (~ 770km s"' ). In this scenario, NGC 5850 would have 
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Fig. 5. — (upper) Temperature map constructed using the mean energy in 
the 0.7— 1 .2 keV energy band (Fe peak) calibrated to temperature using an ab- 
sorbed single APEC model with Galactic absoiption and abundance fixed to 
0.3 Zq . Black lines are contours from the fractional difference surface bright- 
ness image (see Fig. |4), while the dashed ellipse traces the outer (19.8 kpc) 
northeastern and inner (11 kpc) southwestern surface brightness edges shown 
in the lower panel, (lower) 0.75- 2 keV background subtracted, exposure 
corrected merged image of NGC 5846, with the dashed ellipse tracing the 
northeastern and inner southwestern edges. Both the temperature map and 
surface brightness maps are matched in WCS coordinates. Note that in both 
panels, the nucleus of NGC 5 846 is located east of the circular region de- 
noting the ghost bubble, and is significantly offset to the southwest from the 
center of the ellipse. 



approached in front of NGC 5846 from the west, measured 
relative to the elliptical galaxy, and is now receding behind 
NGC 5846 to the east. 

3.2. Gas Temperatures 

We need to understand the temperature structure of the gas 
to understand the thermodynamic state of the gas and deter- 
mine relative gas pressures and velocities. We first extract 
a mean spectrum for NGC 5846 from OBSID 7923 using a 
168" circular region centered at NGC5846's nucleus (RA, 
DEC = 15''06"'29.284^ +0r36'20"25) and local background 
from a circular annulus also centered at the galaxy nucleus, 
with inner, outer radii of 300", 378", respectively. We model 
the spectrum using a VAPEC model with fixed Galactic ab- 
sorption (4.24 X 10^" cm"^), and allow O, Mg, Si, and Fe 
abundances to vary independently. All other abundances are 
fixed at a common fixed value A = 0.5Zq. We fit the spectrum 
over the 0.5-2.0keV energy range to optimize the signal rel- 
ative to background. Although we find a mean temperature 
of 0.65 ± 0.01 keV, and O, Mg, Si,andFe abundances of 0.01, 



0.34, 0.45, and O.26Z0, respectively, this single temperature 
VAPEC model does not provide a statistically acceptable fit to 
the data (x^/dof = 218/95), suggesting the presence of multi- 
temperature gas and strong temperature gradients. 

3.2.1. Temperature Maps 

We investigate these temperature gradients by constructing 
a temperature map for the region. From the azimuthally aver- 
aged global temperature, we know that the temperature of the 
ISM and group IGM gas are both cool. The thermal spectrum 
for gas temperatures < 1 keV is dominated by line emission 
from the the FeL complex. For the ACIS detector, there is a 
nearly linear correspondence between the mean energy mea- 
sured in the FeL peak (0.7 ^ E < 1.2keV) and the temper- 
ature of the gas, as measured by a single temperature APEC 
model (see, e.g. David et al. 2009; Machacek et al. 2010, Fig. 
9). We exploit this correspondence to construct a temperature 
map for the core of the NGC 5846 group. We first calibrate 
the temperature of the gas with the observed mean energy in 
the 0.7- 1.2keV bandpass by fitting a series of simulated ab- 
sorbed APEC model spectra, generated from the spectral re- 
sponse files for our observation, for which the mean energy 
in the Fe peak increases from 0.3 keV to 2keV in 0.05 keV 
intervals. We assume Galactic absorption and fix the abun- 
dance at an average value of O.3Z0. Although this is higher 
than the ^Q.IZq metallicities observed near r^oo in the intra- 
group gas in other cool groups (Rasmussen & Ponman 2007), 
and lower than the mean metallicity ^ 0.5 Zq for ISM gas in 
elliptical galaxies, the temperature is insensitive to changes in 
the abundance over this range, so that this does not signifi- 
cantly affect our results. We average the energy over regions 
determined by adaptively smoothing the 0.7- 1.2keV Chan- 
dra image of the NGC 5846 group and use the derived cali- 
bration curve to associate the mean energy in each region to 
the corresponding temperature of the diffuse gas. Our results 
are shown in Figure |5] 

To check these results and investigate possible tempera- 
ture structure within the central region of NGC 5846, we con- 
struct a second temperature map using the methods of Ran- 
dall et al. (2009) to fit a single temperature APEC model with 
Galactic absorption to regions, containing a minimum of 1000 
counts, grown around each pixel. The abundance is allowed to 
freely vary. The results are shown in the top panel of Figure|6l 
The temperature maps in Figures|5]and|6]are remarkably simi- 
lar, suggesting that the features identified are robust. To guide 
the eye in both maps, we superpose the contours from the 
fractional difference map that shows the spiral structure in the 
X-ray surface brightness (Fig|4]l. The dashed ellipse in both 
maps traces the outer northeast and inner southwest surface 
brightness edges. Note that the nucleus of the galaxy, corre- 
sponding to the X-ray bright AGN, is offset from the center of 
the ellipse. Lower temperatures are found interior to this el- 
lipse and to the northeast of the nucleus, confirming that these 
edges are cold fronts. The gas temperature is again higher 
(^ 1 keV) on the outside of the outer southwestern edge, con- 
firming that this edge is also a cold front. Most of the lowest 
temperature 0.6 keV) gas resides northeast of the nucleus, 
along the contours of maximal fractional difference from the 
best fit /3-model for the X-ray surface brightness. A finger 
of slightly higher (~ 0.7-0.8keV) temperature gas is drawn 
into the nuclear region. Cluster simulations also show that 
sloshing bulk gas motions create multiple edges and may push 
higher temperature gas into the cooler cluster cores. Note, 
however, that there is no evidence for spiral structure in these 
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best-fit APEC model normalization K, where 
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Fig. 6. — Temperature (upper), and pseudo-entropy (lower) maps of 
NGC 5846 constructed using spectral regions, constrained to contain a min- 
imum of 1000 net (background subtracted) counts and fit with single tem- 
perature APEC models using fixed Galactic absorption (4.24 X 10"" cm"^ ). 
Black contours ai'e from the fractional difference surface brightness image 
(Fig. 13), while the white dashed ellipse traces the bright northeastern and 
inner southwestern surface brightness edges, as in Fig. [5| 



temperature maps corresponding to the large scale spiral-like 
northwest extension seen in the X-ray surface brightness. This 
is in sharp contrast to the cool, spiral-like tails observed in the 
NGC 5098 (Randall et al. 2009) and NGC 6868 (Machacek 
et al. 2010) groups, and would be difficult to explain if the 
sloshing encounter occurred in the plane of the sky (Ascasibar 
& Markevitch 2006; ZuHone et al. 2010), again suggesting 
projection effects are important. 

A map of 'pseudo' entropy {s oc TiT^I^) is shown in the 
lower panel of Figure |6l This map is constructed by using the 



as a proxy for for each of the spectral regions in the temper- 
ature map. Da and z are the angular size distance and redshift 
of NGC 5846, respectively, and for simplicity we assume a 
constant emission volume V for each region. As expected, 
the lowest entropy gas is found coincident with NGC5846's 
central region. Entropy increases smoothly with radius. How- 
ever, the distribution is not azimuthally symmetric, but again 
follows the off-centered ellipse (dashed line) that traces the 
northeastern outer and southwestern inner edges. While this 
map is suggestive, it provides only a qualitative picture for 
entropy throughout the region. First, the uncertainties in the 
abundances for spectral regions containing only 1000 source 
counts are large. Although the model temperatures do not 
vary significantly with abundance in this range, the normal- 
ization K in equation [T] does by a factor of > 2 for 0.64 keV 
gas and abundances between 0.4- \ Zq (see, e.g. Machacek 
et al. 2006). This increases the uncertainty in the derived den- 
sity to ~ 40% and in the derived entropy to ~ 26%. Sec- 
ond, the true geometries (volumes) of the emission regions 
in the disturbed gas are likely complex, as demonstrated by 
gas distributions in the cores of simulated clusters (see, e.g. 
Ascasibar & Markevitch 2006; ZuHone & Markevitch 2009; 
ZuHone et al. 2010), and may not be well represented by the 
simple model used to construct the maps. Thus the qualita- 
tive features shown in these maps need to be checked by more 
targeted spectral analysis and modeling. 

3.2.2. Modeling Abundances 

Since we want to constrain the relative densities between 
differing spectral regions as well as possible, we choose to 
use the VAPEC thermal plasma model for detailed spectral 
fits. To optimize signal to noise, the fit range is restricted 
to the 0.5-2keV band, where emission from cool < Ikev 
gas dominates. Guided by the temperature maps in Figures 
|5] and |6] we extract the spectrum for an X-ray luminous re- 
gion expected to be homogeneous in temperature, defined 
as the sum of two circular regions with radius of 39" cen- 
tered at (15''06"'29.85', +01°37'14f'89) and 31."4 centered 
at (15''06'"34.445'', +01°37'38f'50), respectively. We model 
the spectrum of this region using an absorbed VAPEC model 
with fixed Galactic absorption, but with the abundances for 
O, Mg, Si, and Fe, expected to be important for low tem- 
perature gas, free to vary. The remaining abundances for 
elements that do not have strong emission features in the 
0.5-2keV energy band are fixed at Q.5Zq. We find a temper- 
ature of 0.60 ± 0.01 keV and O, Mg, Si, and Fe abundances 
of O.24![!;|lZ0, O.53t!^Z0, Q.ll^f^ZQ, and O.49![!;i^Z0, 
respectively (x^/dof = 79.8/79). We then fix the metal abun- 
dances at these values in our spectral models for the less lu- 
minous features. 

3.2.3. Azimuthally Averaged Temperature Profile 

In Figure |2l we show the azimuthally averaged temper- 
ature profile of NGC 5846 to a maximum radius of 40kpc 
using 14 circular annuli, with logarithmically increasing bin 
width ranging from 6" (0.7 kpc) for the innermost bin to 64 ."2 
(7.45 kpc) for the outermost bin (see Table[T]). The local back- 
ground for all bins was chosen to be a 3' radius circular region. 



Edges and Bubbles in NGC 5846 



7 





1 






1 ^ 

I - 












\ i 





















Fig. 7. — Projected, radial temperature profile of the NGC 5846 group. 
The spectra for logarithmically increasing bin widths are modeled using an 
absorbed VAPEC plasma model with fixed Galactic absorption (see Table 
[T). Vertical lines denote the positions of e dges to the northeast (dashed) and 
southwest (dotted), (see Figs.|2]and|5]and j|3.3.U . 



centered at RA, Dec of 15''05'"52.374', +0r35'43."9, placing 
the background region as far as possible (9.'2) from NGC 5846 
to minimize the contribution from group emission, while still 
being completely contained on the ACIS-I detector. 

We model the spectra with a single absorbed V APEC 
model, Galactic absorption and abundances fixed as in ^3.2.21 
and fit the spectrum over the 0.5-2keV range. Our results are 
listed in Table [T] A spectral fit is deemed statistically accept- 
able if the probability of getting a as large or larger than ob- 
served if the model is correct, denoted the 'null' parameter, is 
null > 0.05. Although the single temperature VAPEC model 
provides a reasonable description of the spectra outside the 
central 1 .5 kpc of the galaxy, the x^/dof is only marginally ac- 
ceptable (null ^ 0.01 -0.02) for bins with mean radii of 8 kpc, 
16 kpc and 25 kpc, and the residuals in those bins suggest 
the presence of multi-temperature gas (Buote 2000). This is 
not surprising given the elliptical, off-centered distribution of 
temperatures shown in the temperature maps. Within lOkpc, 
the gas temperature is roughly isothermal with temperatures 
of 0.62- 0.65 keV. The mean temperature rises modestly to ~ 
0.75 keV at ^ 10 kpc and then rises steeply to > 1 keV at radii 
beyond 20 kpc, consistent with ROSAT and XMM-Newton re- 
sults (Trinchieri et al. 1997; Finoguenov et al. 1999; Nagino 
& Matsushita 2009). The temperature increases for increas- 
ing radii at 8 kpc and 20 kpc are approximately coincident 
with the positions of the inner and outer edges seen in the 
surface brightness profiles (Fig. |8]l. In the central 0.7 kpc, a 
single temperature VAPEC model, with best fit temperature 
0.72 ±0.02 keV (xVdof= 76/48, null = 0.006) is not a good 
description of the data, with the see-saw pattern in the residu- 
als indicating the presence of multi-phase gas. A t wo tem per- 
ature VAPEC model with abundances fixed as in 33.2.21 and 



temperatures 0-60^'^, keV and 1 .3!!!;;^ keV (xVdof = 3 1 /46, 
null = 0.95)is a much better description of the data. The single 
temperature VAPEC model also fails to describe the spectrum 
in the 0.7- 1.4 kpc annulus. This central region r < 1.5 kpc) 
is complex, showing X-ray cavities (bubbles) filled with ra- 
dio plasma from recent AGN activity, and X-ray bright knots 
as well as Ha emission (Fig. [12] and ® . The presence of 
higher temperature gas in the spectral model for the central 
region of NGC 5846 is unusual for group dominant elUptical 



galaxies, whose bright cores tend to be cool. This spectral 
component may be an artifact of applying a simple spectral 
model to multi-phase gas in the central region of NGC 5846, 
or indicate the presence of higher temperature gas that has 
been heated by the recent passage of a shock. 

3.3. Modeling Gas Properties Across the Surface 
Brightness Edges 

3.3.1. Gas Densities 

To study the northeastern edges quantitatively, we identify a 
bounding ellipse centered at NGC 5846's nucleus with semi- 
major (-minor) axes of 2'A9 (2fl6) and position angle 42.3°, 
that traces the outer sharp surface brightness discontinuity in 
the angular sector extending from 123° to 180°. We con- 
struct the surface brightness profile, in that angular sector, 
using concentric elliptical annuli of increasing (decreasing) 
logarithmic width for radii outside (inside) the bounding el- 
lipse. To the southwest, we construct the radial surface bright- 
ness profile using concentric circular annuli, centered again 
on NGC 5846 and with logarithmically increasing width, con- 
strained to lie in the angular sector extending from 268° to 
338°. In the right panel of Figure |3] we compare the mean 
/3-model (dashed line) with these surface brightness profiles, 
confirming the qualitative evidence for multiple edges seen in 
Figure 12] Not only is the surface brightness and thus gas den- 
sity much higher at < 20 kpc inside the outer northeast edge, 
as also seen in ROSAT HRJ images (Trinchieri et al. 1997), 
but the edge is sharper and deeper than that to the southwest, 
suggesting that the ISM gas distribution has been displaced 
significantly to the northeast. The surface brightness profile 
to the southwest is a factor 2 fainter at ^ 10 kpc compared to 
that at the same distance from the nucleus to the northeast. 
There is a shallower discontinuity near 20 kpc in the south- 
western profile than in the northeast, perhaps a remnant from 
a previous oscillation. 

Following Vikhlinin et al. (2001), we denote the edge posi- 
tion and ratio of electron density across the inner (outer) edge 
in each sector as r\ and J\ (r2 and J2 ), respectively. The elec- 
tron density on either side of each edge is modeled with a 
power law. Thus, 



ne = J\ 



«o/2(;;^)""](^))"'(r<n) (2) 

/ \ amid 

= «o/2 7: {ri<r<r2) (3) 



:«() 



(r > r2) 



(4) 



where the power law slopes of the density for radii inside the 
inner edge, between the edges, and outside the outer edge are 
Q^b c^mid^ and a2, respectively and «o is a constant normaliza- 
tion. We then integrate the electron density model along the 
line of sight and use a multivariate minimization scheme 
iteratively to fit the model parameters for each sector We first 
determine a2 and r2 from a fit to the outer edge alone, and then 
fix these parameters in a fit across both edges while allowing 
the remaining parameters (7i, J2, r\, ai, amid, no) to vary. Our 
results for the northeast (NE) and southwest (SW) sectors are 
listed in Table [2] and shown as solid curved lines in the up- 
per and lower panels of Figure [8] Edge locations are denoted 
by vertical lines. In the northeastern sector, we find a sharp, 
deep outer edge located 19.8 kpc from the nucleus with a gas 
density jump, assuming no strong abundance changes at the 
boundary, of 2.9!j]:i The inner edge, located at 6.6 -6.8 kpc. 
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Fig. 8. — (upper) 0.5 — 2keV surface brightness profile to the northeast of 
NGC 5846, taken in elliptical annuli concentric to the bounding ellipse, cen- 
tered at NGCC5846's nucleus, with semi-major (-minor) axes of 2'A9 (2f 16) 
and position angle 42.3° and constrained to lie in the angular sector between 
123° to 180°. The solid line represents the best-fit density law model from 
equations 1 -|4] and listed in Table |2] The inner (outer) edge location at 
~ 6.7 kpc (19.8 kpc) is denoted by the dotted (long dashed) vertical line, re- 
spectively, (lower) 0.5 — 2keV surface brightness profile to the southwest 
of NGC 5846, taken in concentric circular annuli constrained to lie within 
the angular sector between 268° to 338°. The solid line is the best-fit den- 
sity model given in Table |2] The inner (outer) edge locations at 11.0 kpc 
(19.1 kpc) are denoted by dot-dashed (dashed) vertical lines, respectively. 
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Fig. 9. — (upper) Projected, radial temperature profile across the northeast 
surface brightness edges. Spectral regions are elliptical annuli chosen to be 
concentric to the bounding ellipse defined in Figjs] and constrained to lie in 
the angular sector from 123° to 180° (see Table|3)7 (lower) Projected, radial 
temperature profile across the southwestern surface brightness edges. Spec- 
tral regions are concentric circular annuli centered on NGC 5846's nucleus 
and constrained to lie within the angular sector from 268.2° to 338° (See 
Table|4). For each case vertical lines denote the edge positions, and the spec- 
tra are modeled using an absorbe d VAPE C plasma model with fixed Galactic 
absorption and abundances as in j|3.2.2l and Table[T] 



is not as well-defined and is a factor ^ 2 shallower with a 
density jump of 1.4 ±0.2. Similarly we find edges at 19.1 kpc 
and 11.0 kpc from the nucleus in the southwest sector with 
density discontinuities of 2.1;I;q 4 and 1.6 ±0.2, respectively. 
The inner edge in the southwest sector falls close to a chip 
gap on the ACIS-I detector. To test whether this might bias 
our results, we fit the surface brightness profile from the same 
SW sector using only data from the 30 ks obsid 788 observa- 
tion, where the inner edge region of interest falls entirely on 
the S3 CCD, and found no significant differences in the fit re- 
sults. We also note that the /3-model parameters derived from 
the outer power law slopes (02) for the NE sector (0.39;!;y 
and SW sector (0.39;!;Qog) are in good agreement with each 
other and with the /3 = 0.45 obtained from a single /3-model 
fit to the NGC 5846 group's mean (azimuthally-averaged) ra- 
dial surface brightness profile, shown in the upper panel of 
Figure |3] 

3.3.2. Gas Temperatures, Pressures and Velocities 



To better understand the temperature changes across the 
surface brightness edges, we extract spectra from either side 
of the edge using regions concentric with the edge morphol- 
ogy. For the northeastern edges, we extract spectra from el- 
liptical annular regions concentric to the same bounding el- 
lipse that was used to generate the northeastern sector surface 
brightness profiles and constrained to lie in the same angular 
sector (see ED and Fig. |8]l. Regions NEO-1 (NEOO) fall 
just inside (outside) the outer northeastern edge at 19.8 kpc. 
NEI-1 (NEIO) lie inside (outside) the inner northeastern edge 
at ~ 6.7 kpc. Similarly for the southwestern edges, we use 
concentric circular annuli, centered on NGC5846's nucleus 
and constrained to lie in the same sector used to construct the 
southwestern surface brightness profile, to measure the gas 
temperatures on either side of each edge. Regions SWO- 
l(SWOO) fall just inside (outside) the outer southwestern 
edge at 19.1 kpc, and SWI-1 (SWIO) lie just inside (outside) 
the inner southwestern edge at 11.0 kpc, respectively. Back- 
grounds are taken from renormalized blank sky backgrounds 
provided by the Chandra X-ray Center for the same regions, 
normaUzed to the observations by comparing count rates in 
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Fig. 10. — 610 MHz radio contours superposed on the background sub- 
tracted, exposure connected, coadded 0.75 — 2 keV Chandra image of diffuse 
gas in NGC 5846. Anows show the location of the inner northeastern edge. 



the 9.0- 12.0keV energy band (see ©. 

Our results are listed in Table [3] for the northeastern edges 
and in Table |4] for the southwestern edges, and are plotted in 
Figure |9] As expected, the temperature profiles are not sym- 
metric. Lower temperature gas has been pushed to larger radii 
to the northeast. For the outer northeastern edge, the temper- 
ature rises from A:7} = 0.7 1 ± 0.02 keV inside the edge to kTo = 
1.00 ±0.04 outside the edge, such that the temperature jump 
between undisturbed gas outside the cold front and gas just 
within the cold front edge is Ti/To = 0.71 ±0.05. Given the 
density jump of n,/no = 2.9;!;y 4 from Table|2] we find a pres- 
sure ratio Pi/ po = 2. 1^0 4 across the edge. If this pressure ratio 
is interpreted as evidence for bulk motion of the cold front, 
then, following Vikhlinin et al. (2001), we infer that the front 
is moving transonically with velocity ^ 520;!;jo()kms"' (Mach 
1.0 ±0.2) relative to the 1 keV gas outside the cold front. 
Thermal pressures across the inner edge are, within the large 
uncertainties, consistent with pressure equilibrium and zero 
relative velocity. Although the uncertainties are larger, the 
pressure ratios derived across the southwestern inner (outer) 
edges are remarkably similar to those to the northeast. The 
pressure ratio across the inner southwestern edge is consis- 
tent with thermal pressure equilibrium, while the pressure ra- 
tio across the outer southwestern edge is ~ 2, suggesting pos- 
sible transonic motion of the gas. 

In the upper panel of Figure |9] we also see that there is 
no statistically significant change in the temperature across 
the northeastern inner edge. Since these are projected tem- 
peratures, higher temperature gas at large radii along the line 
of sight may reduce the observed projected temperature ratio 
over that of the true (deprojected) ratio of temperatures across 
the front. 

However, the origin of the inner northeastern edge may 
be more complex than simple gas sloshing. In the radio 
band, NGC5846 hosts a low radio power (1.5 x 10^' WHz"' 
at 1.45 GHz) source (see also Fig. [T2] |. Giacintucci et al. 
(2011) observed the source with the Giant Metrewave Radio 
Telescope at 610MHz. Their 15"-resolution image detected 
a central unresolved component and an extended feature of 
only ~ 3 kpc toward the northeast. In Figure [TO] we super- 
pose the 6 10 MHz radio emission contours from Giacintucci 
et al. (2011) on the 0.75-2.0kev image of NGC5846, and 
find that the inner northeastern edge lies, in projection, just 



North Bubble 




Fig . 1 1 .— 0.5 - 2 keV close-up view of the central ~ 10 kpc of NGC 5846 
showing the ghost cavity located ~ 5.2 kpc from the nucleus, and the north 
and south partially disrupted inner X-ray bubbles with centers ~ 0.6 kpc from 
the nucleus. 



outside the extended structure detected at 6 10 MHz. This sug- 
gests the possibility that the inner northeastern edge may be 
cool gas displaced by radio plasma in a bubble formed during 
a prior episode of AGN activity. 

4. X-RAY CAVITIES AND NUCLEAR ACTIVITY 

The thermodynamic state of gas in the central 2 kpc of 
NGC 5846 is complex, with a non-radial, non-azimuthally 
symmetric distribution due both to merging activity and 
episodic feedback from the central AGN. Observations in 
multiple wavebands (radio, optical and X-ray) show that the 
gas is truly multi-phase with a mix of cool, dusty filaments, 
warm gas observed through its optical line emission, and hot 
gas with a range of temperatures and densities. The gas is dis- 
turbed both by merging activity and ongoing AGN outbursts, 
resulting in observed non-radial, non-azimuthally-symmetric 
gas distributions that are more difficult to interpret. 

Evidence for two recent AGN outbursts are seen in the X- 
ray 'ghost' cavity found 5.2 kpc east of the nucleus, shown 
in Figure |2] and in two roughly symmetrical, inner cavities 
(bubbles), observed previously by Trinchieri & Goudfrooij 
(2002) and Allen et al. (2006), whose centers lie 0.6 kpc 
from NGC5846's nucleus (see Figures [TTI and [T2]|. The in- 
ner bubbles are highly significant, with the excess flux in the 
rims more than lOcr greater than that expected from back- 
ground. While the rims on the bubbles do not appear com- 
plete (particularly to the south), it is unclear whether this is 
the result of actual bubble disruption through the onset of hy- 
drodynamic instabilities or a projection effect caused by the 
evolution of the bubble through multi-phase gas, already dis- 
turbed by merging activity. We have retrieved and reanalyzed 
1.45 and 5 GHz data from the VLA archive (see Table |6] for 
details). Calibration and imaging were carried out using the 
NRAO Astronomical Image Processing System (AIPS) using 
standard methods. We find that radio emission at 1.45 GHz 
and 5 GHz is coincident with the inner bubbles, and that the 
1 .45 GHz radio emitting plasma appears to extend outside the 
apparently roughly spherical bubble walls to the northeast and 
south. However, no radio emission is detected at 6 10 MHz 
from the ghost cavity to the west. 

In the upper right panel of Figure [121 we show contours 
from an Hq;-i-[NII] image taken with the 4.1 m Southern Ob- 
servatory for Astrophysical Research telescope. Hq;-i-[NII] 
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Fig. 12. — A multiwavelength .study of the central 4kpc region of NGC5846, (upper left) 0.5-2keV exposure corrected, coadded Chandra X-ray image 
superposed with 1 .45 GHz VLA radio contours. 1 pixel = 0"5 X 0"5 and the image has been smoothed with a 1 " Gaussian kernel. The radio beam has 6" X 5" 
beamsize. Contours begin at 0.1 1 mJy/beam and increase by factors of 2.(upper right) Same Chandra image superposed with higher resolution 1.45 GHz VLA 
radio contours to highlight the small scale radio structure. The radio beam has 1.5" X 1.5" beamsize. {lower left) The same Chandra image superposed with 
5 GHz VLA radio contours. The radio beam has 9" X 5" beamsize. Contours begin at 0.07 mJy/beam and increase by factors of 2. {lower right) The same 
Chandra image supeiposed with Ho -l-[NII] emission contours, from an image taken with the 4.1 m Southern Observatory Astrophysical Research telescope. 
Note the close correspondence with small scale structure in the 1 .45 GHz emission. 



emission is weak, with a total Ha flux of ^ 6.6 x 
10"'"^ erg s"' cm~^ . We assumed the following ratios between 
the [Nil] Unes and the Ha line: [Nll]6583/Ha = 1 and 
[N11]6548/[N11]6583 = 0.35. With the Kennicutt (1998) re- 
lation, where the star formation rate is given by 



SFR(Moyr-') = 7.9 : 



10-''(Lhc 



)ergss 



(5) 



the star formation rate in NGC 5846 is only 0.037 M© yr"' . 
The HaH-[Nll] filaments seem dusty and are aligned with the 
X-ray substructures. Ha emission is found inside the inner 
edge of the northern bubble rim, and also just inside the visi- 
ble rim of the bubble to the south. 

4. 1 . AGN Outburst Energetics 

We estimate the energy of the outburst by the enthalpy 
i^pV) carried in the cavity, assuming that the cavity is filled 
with relativistic (7 = 4/3) plasma, p is the mean thermal 
pressure in the undisturbed ambient gas at the same radius as 
the center of the bubble, and V is the bubble volume (Birzan 
et al. 2004, Allen et al. 2006). The undisturbed gas pressure 
is estimated using the mean /3-model fit to the density from 
ED and temperatures of 0.65!!]:jj^keV and 0.78![i:jj| for the 



inner and ghost bubbles, respectively, taken from the temper- 
ature maps at the bubble locations. For simplicity we assume 
each bubble is spherical. The method used to estimate the 
bubble age depends on the evolutionary stage of the bubble. 
For young bubbles that are still momentum driven by the ra- 
dio jets, the age of the bubble is taken to be the time f, = d/cs 
for the bubble to rise to its given location d, where d is the 
projected distance of the center of the bubble from the galaxy 
nucleus, at the speed of sound Cj in the ambient gas (Omma 
et al. 2004; Birzan et al. 2004). For more evolved bubbles, the 
age of the bubble may be better approximated using the time 
(tb) for the bubble to rise buoyantly at its terminal velocity, 
tb = d(SC/2gVy^^, where S is the bubble cross-section, V is 
the bubble volume, C = 0.75 is the drag coefficient, and g is the 
acceleration of gravity at the projected distance d of the bub- 
ble center from NGC5846's nucleus (Churazov et al. 2001; 
Birzan et al. 2004). Once the bubble enthalpy and age are 
determined, the instantaneous mechanical power of the out- 
burst is given by Lmech = E /t. Our results are given in Table 
|7] where we have used the buoyancy timescale ti, to estimate 

^mech- 

For the inner bubbles, we find outburst energies ^ 
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Fig. 13. — Unshaip masked 0.5-2 keV background-subtracted, exposure- 
corrected, coadded Chandra image of the inner bubbles in NGC 5846. The 
unsharp masking used Gaussian smoothing scales of 1 and 45 pixels with an 
exposure map threshold of 20%. 1 pixel = Of'5 X 0'.'5. 

10^^ ergs, bubble ages ^ 2Myr and mechanical powers ~ 
2-4 X 10'*' ergs"', consistent with previous work (Allen 
et al. 2006). The outburst energy inferred for the 'ghost' bub- 
ble is modestly larger 5 x 10^^ ergs) and the bubble age is 
~ 12Myr, giving the time between outbursts (AGN duty cy- 
cle) ~ lOMyr. However, the uncertainties in these measure- 
ments are large. First the three dimensional cavity geometry is 
not known and may affect estimates of the bubble volume and, 
consequently, outburst energy by factors of a few. Second, our 
simplified model assumes AGN bubbles rising in a static at- 
mosphere. However, both the inner and ghost bubble ages are 
shorter than the expected gas sloshing timescales given by the 
time since closest approach of the proposed galaxy perturber 
NGC 5850 (60-200Myr; Higdon et al. 1998). Thus non- 
hydrostatic bulk gas motions may have influenced AGN ac- 
tivity and likely also affect bubble morphology, the observed 
distance of the bubble from the nucleus, and the probability 
for bubble disruption. Numerical simulations are required to 
model bubble evolution in this complex, dynamical environ- 
ment. 

4.2. Knots 

In Figure [T3] we use unsharp masking to reveal the inter- 
nal structure of the northern and southern inner bubbles in 
NGC 5846. The unsharp masked X-ray image is created using 
Gaussian smoothing scales of 1 and 45 pixels, respectively, on 
the full resolution exposure corrected, co-added 0.5-2keV 
Chandra image. An exposure map threshold of 0.20 was im- 
posed to maximize signal-to-noise in our result. We see that 
the inner bubble rims are not smooth, but are threaded by 9 
knot-like features (see Table [8j. To determine the detection 
significance of these knots, we calculate the cumulative prob- 
abihty that the observed counts or higher in each knot in the 
0.5— 2keV energy band could have arisen from aPoisson fluc- 
tuation of the surrounding rim emission. Knots 4 and 8 are 
least significant with a cumulative probability of 0.02, while 



for all other knot features, the cumulative probability that the 
observed counts or higher might be a fluctuation of the sur- 
rounding rim emission is < 1.7 x 10"'*. Thus the significance 
of the knots as distinct features is high. 

We want to compare the spectra of these knot-like struc- 
tures to the spectrum of diffuse gas in the surrounding rims. 
In Table |8] we define circular spectral regions around each 
knot and elliptical annuli excluding the knots and nucleus that 
trace the north and south bubble rims. An elliptical annu- 
lus (excluding point sources) located just outside the rims is 
taken as the local background for both the knot and rim re- 
gions. We model the spectrum of the rims with an absorbed 
VAPEC model ass uming Galactic absorption and fixing the 
abundances as in ^3.2.21 As shown in Table [TT] the mean 
temperature of gas in the north rim is QJl^QlksY for the 
89.5 ks exposure (OBSID 7923). However, the model does 
not provide a statistically acceptable fit (x^/dof = 111/65, 
null = 0.03), and the residuals display the characteristic see- 
saw pattern of multi-temperature gas (Buote 2000). In con- 
trast, hot gas in the south rim region is well modeled by an 
absorbed single temperature VAPEC model. We find a best fit 
gas temperature of 0.64+|j{]3 keV (x^/dof = 94.3/97), where 
we have used a simultaneous fit to both the 89.5 ks (OBSID 
7923) and 23.2 ks (OBSID 788) data sets to reduce uncertain- 
ties. The more complex spectral structure for gas in the north 
rim compared to that in the south rim region may be due to 
projection effects, such that multi-phase gas is more promi- 
nent along our line of sight to the north than to the south. 
Alternatively, the combination of AGN and merging activity 
may have introduced hydrodynamic instabilities in the fluid 
that were more efficient at disrupting the bubble to the south 
than to the north. The intrinsic 0.5-2keV luminosity of dif- 
fuse gas in the north rim (1.0 x lO'^'ergs"' ) is 50% higher 
than that in the south rim (0.65 x 10""' ergs"' ). 

The individual knot-like features are X-ray faint. In Ta- 
ble |9] we list the net source counts for the individual knots in 
four energy bands: 0.6 - 0.8 keV (FeSS), 0.8 - 1 . 1 ke V (FeS), 
1.1 -2.0keV (M), 2-4keV (H). The first two bands ai-e cho- 
sen to divide the energies containing the FeL line complex 
(the Fe peak from 0.6- 1.1 keV) in half, such that the hard- 
ness ratio constructed from these bands will be sensitive to 
shifts in the energy of the peak emission from the FeL com- 
plex, caused by temperature differences in cool (< 1 keV) gas. 
For comparison, we also include net source counts in these 
energy bands for the north and south rims. Notice that in all 
cases (knots and rims), the X-ray emission is dominated by 
the Fe peak, confirming that the X-ray emitting gas is cool, 
and, for all the knot-like features except Knot9, the net photon 
counts above 2 keV are consistent with zero at the < 2cr con- 
fidence level. Since the number of photons in each individual 
knot is small, we use the hardness ratio across the iron peak, 
(FeS -FeSS)/(FeS + FeSS), to search for temperature differ- 
ences between the individual knots and between the knots and 
the rims. We also compute the standard medium band hard- 
ness ratio (M - S)/(S + M+H), where S is the number of pho- 
tons in the 0.6- 1.1 keV energy band. These hardness ratios 
along with the total (0.6-4keV) net source counts for each 
region are given in Table [TOl Although the hardness ratios for 
Knot4 and Knot5 suggest a mildly softer spectrum than the 
others, the most striking result from Table [Tolis that the hard- 
ness ratios of the knots are nearly identical to that of the rims. 
We stack the knot spectra to determine quantitatively a mean 
temperature for the knots, and confirm that the best fit temper- 
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ature {0.12]^^^, x^/dof = 67.3/45 for the 89.5 ks exposure) is 
statistically indistinguishable from that measured in the sur- 
rounding rims. The total 0.5-2keV luminosity of the knots 
is 0.3 X lO'**^ erg s~' , one third of the X-ray luminosity of dif- 
fuse gas in the north rim. 

To investigate the gas density and thermal pressure in the 
knots, we fit the spectrum of the brightest knot (Knot6 in the 
north rim with 0.5-2keV luminosity of 5.7 x lO-'^ergs"' ) 
as an example. We fix the VAPEC model temperature to 
0.72 keV, the best fit temperature for the stacked knot spec- 
trum (and gas in the north rim), and allow only the VAPEC 
model normalization to vary. We then compute the electron 
density from the VAPEC normalization (see eq. [U, assum- 
ing that the knot is a sphere with radius 1"4 (0.162kpc) 
uniformly filled with gas. We find an electron density of 
0.33 ± 0.03 cm"-' and corresponding thermal pressure of 7.3 x 
10"'" ergcm"^ where the uncertainties reflect only the statisti- 
cal uncertainties of the spectral model fit. For such dense gas, 
the cooling time is short, ^ 30 Myr. 

For comparison, we compute the electron density and cor- 
responding thermal pressure for the diffuse gas in the north 
rim surrounding Knot6. Again the largest uncertainty in the 
analysis is the assumed morphology of the emission region. 
To bound the electron density in the rim, we compute the 
electron density assuming first that the observed emission is 
gas uniformly distributed throughout a cylindrical shell with a 
cross sectional area that of the spectral extraction region and a 
height equal to the bubble radius (an upper limit on the emis- 
sion volume and consequently lower limit on the derived den- 
sity) and, second, that the X-ray emission comes solely from 
a spherical shell with dimensions determined by the spectral 
extraction region (a lower limit on the emission volume), ex- 
cluding the knot regions in each case. The electron density in 
the north rim is then bounded by 0.076 < < 0.096 cm"^ . 
This is somewhat larger, as expected by the greater surface 
brightness of the rim, but comparable to the electron gas 
density (0.074 c m"'' ) derived from the mean /3-model for 
NGC 5846 (see ^3.11) . We infer a thermal pressure in the gas 
surrounding Knot6 of (1.7-2) X 10"^°ergcm"^ . atleastafac- 
tor 3 lower than that derived for Knot6, assuming a spherical 
geometry. 

There are two possible scenarios to explain this discrep- 
ancy. If the knots are stable structures in pressure equilibrium 
with their surroundings, then either the knots are not spheri- 
cal, cooling blobs or they are stabilized by non-thermal pres- 
sures, or both. In this scenario, the most likely explanation is 
that the true 3-dimensional morphology of the knot-like struc- 
tures is more extended and complex. For example if they are 
filamentary rather than spherical, i.e. more like the X-ray fila- 
ment observed in M87 (Forman et al. 2007) but viewed with 
the filament oriented along the line of sight, the emission vol- 
ume would be larger, implying a lower electron density and 
allowing the structure to be closer to or in thermal equilib- 
rium with its surroundings. Then the knot-like appearance of 
the structures would be primarily a projection effect. It may 
be difficult to explain, however, why all the knots would have 
this orientation. 

The second possibiUty is that the pressure estimates are, 
in fact, roughly correct and the knots and surrounding fea- 
tures are the visible result of a strong shock being driven into 
the ambient gas by the supersonic inflation of the radio lobe. 
The general morphology of the rim and knots is similar to the 
high Mach number shocks seen in NGC 4552 (Machacek et al 



2006) and Her A (Nulsen et al. 2005). There is a noticeably 
bright, roughly elliptical, enhancement associated with the ra- 
dio lobes in both cases. The morphology of the knots and the 
rim looks similar, albeit at much lower signal to noise. The 
knotty morphology could be the result of a shock being driven 
into a multi-phase ISM. In this scenario, the knots are the re- 
sult of cooler clouds embedded in a hotter medium. A shock 
driven through such a stratified medium will drive spherical 
shocks into the clouds. They are brighter because they are 
denser, while the post-shock gas must be in rough pressure 
equilibrium. 

We can use the density between the rim and the pre-shocked 
gas in the Rankine-Hugoniot conditions to estimate the Mach 
number of the shock and the expected temperature increase 
in the rim (Landau & Lifschitz 1959). Assuming the bub- 
ble rims are spherical shells and th at th e preshock density is 
given by our mean /3-model from ^3.11 evaluated at the mid- 
point of the rim, we can use the density jump between the rim 
and the pre-shocked gas in the Rankine-Huygenot conditions 
to estimate a shock Mach number. From the more complete 
northern rim, we find a density jump of 1.4, and Mach num- 
ber 1.3. However, the Rankine-Huygenot conditions assume 
the shock front is thin. Our measured value for the density 
jump is an average value across the rim, and so will under- 
estimate the true density jump at the shock front. Thus this 
inferred Mach number should be interpreted as a lower bound 
on the true Mach number for the shock. These averaged val- 
ues for the north rim in NGC 5846 are similar to those mea- 
sured for bubbles in the central region of NGC 4552, where 
more detailed modeling indicated a shock Mach number of 
^ 1 .7 (Machacek et al. 2006). The measured (projected) tem- 
perature jump across the northern rim in NGC 5846 (^ 1 . 1) is 
expected to be smaller than the Rankine-Hugoniot prediction 
(1.3) due to dilution by emission from unshocked gas along 
the line of sight. Note that knotty features are also seen in 
X-ray images of NGC4552's shocked bubble rims (see, e.g, 
the right panel of Fig. 1 in Machacek et al. 2006). 

5. NGC5846A: GALAXY INFALL & GAS STRIPPING 

NGC5846A, shown in Figure [U is directly south of the 
dominant group eUiptical galaxy NGC 5846. NGC5846A 
was among the earliest galaxies confirmed to be a compact 
elliptical galaxy (cE). Compact elliptical galaxies form a rare 
class of early type galaxies that share many of the unusual 
properties of M32, companion to M3 1 (Andromeda) in the 
Local Group (Faber 1973). Compact elliptical (cE) galax- 
ies have low optical luminosities {^^ lO^L©), similar to the 
more common dwarf elliptical galaxies, but their half-light 
radii (~' 0.25 kpc) are several times smaller, such that their 
optical surface brightness is much higher than expected for 
their luminosity (e.g. Faber 1973; Smith Castelli et al. 2008; 
Chiligarian et al. 2009). All are found near massive com- 
panions. They also tend to have higher stellar velocity dis- 
persions and higher metallicities than expected for their lu- 
minosities, and lie above the locus of normal dwarf ellip- 
tical galaxies on the Faber-Jackson relation (Faber & Jack- 
son 1976). Currently only about a dozen cE galaxies have 
been found within 100 Mpc, where their properties can be 
confirmed by ground-based telescopes (for recent searches 
see Smith-CastelH et al. 2008, 2010 in the Antha cluster; 
Price et al. 2009 in Coma; Chiligarian & Bergond 2010 in 
the NGC 5846 group). 

While the origin of cE galaxies and their relationship to 
galaxy evolution remain under debate, the most likely expla- 
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Fig. 14. — (upper) 0.3 — 2k:eV background-subtracted, exposure-coiTected, 
coadded Chandra image of the center of NGC 5846 showing the infalling cE 
galaxy NGC 5846A. ( lower) Close-up of the image in the left panel highlight- 
ing NGC 5846A's small gaseous corona and stripped tail. 1 pixel = 0"5 X Q'.'5 
and the image has been smoothed with a 1"0 Gaussian kernel. 



nation for their high surface brightness, high velocity disper- 
sion and compact structure is that their progenitors were mas- 
sive, possibly disky, galaxies that have undergone severe tidal 
threshing from interactions with a nearby massive compan- 
ion. The outer stellar layers have been tidally stripped from 
these galaxies leaving only the compact stellar core (Bekki 
et al. 2001, 2003). The near-infrared indices for NGC5846A 
suggest that the progenitor of NGC5846A may have been a 
massive elliptical galaxy, and that the catastrophic tidal in- 
teraction, responsible for stripping ~ 90% of its progenitor's 
stellar mass, may also have funneled gas into NGC5846A's 
core, resulting in an intermediate age population of stars near 
its nucleus (Davidge et al. 2008). 

One might expect such strong tidal interactions would 
also remove any hot gas corona associated with the com- 
pact elliptical galaxy, leaving it devoid of X-ray emitting 
gas. Figure [14] shows, however, an X-ray tail extending 
0.54 kpc (4 ."7) to the south of NGC5846A, evidence for on- 
going ram pressure stripping of a residual gaseous corona 
of the compact elliptical galaxy. Since the radial veloc- 
ity difference between NGC5846A and NGC 5846 is - 
487 km s"' (NED), already transonic for ~ 0.7 keV gas, the 
infall is likely supersonic. Assuming a thermal model with 
kT ^ 0.5 - 1 keV, uniform filling and cylindrical geometry 
for the tail, the mean density of gas in the tail is estimated 
to be ^ 0.05 cm"-' and mass in the tail of ~ 10^ Mq. The 
average mass return rate for a simple population of stars in 
a galaxy of age f», where f, is the end of the galaxy for- 
mation epoch and stars are passively evolving, is dm^/dt ~ 



1.5 X lO-"(r*/15)"'-^LBM0yr-i (Ciotti etal. 1991). Assum- 
ing that NGC5846A shares many of the same properties as the 
compact elUptical galaxy NGC4486B (Davidge et al. 2008), 
such that ^ 12 Gyr (Soria et al. 2006), and has an absolute 
B-band magnitude ~ -17 (Davidge 1991), the rate of mass 
loss by aging stars is ^ O.O18M0yr~'. The time between 
stripping episodes needed to replenish the ISM with sufficient 
mass to explain the gas mass observed in the tail is '--^ 5.5 Myr. 

The central X-ray corona in NGC5846A has a radius ^ 
1"3, also of the order of the Chandra point spread func- 
tion. Since this is consistent with a point source, we can 
only set an upper limit of < 150pc for the physical radius of 
NGC5846A's gas corona. Interpretation of the X-ray prop- 
erties of the central region of NGC5846A is further com- 
plicated by the presence of a nuclear point source. The ob- 
served 2- lOkeV source count rate in a 1"3 circular region 
surrounding the nucleus (~ (3.4 ±0.6) x lO"'* counts s"' ) is 
more than an order of magnitude larger than that expected 
if X-ray emission in the region was due solely to hot the X- 
ray emission from the region were due to ~ 1 keV hot gas 
with solar abundance. Assuming a power law model with 
photon index F = 1 .4 for the nonthermal emission, we find 
a 2- 10 keV nuclear luminosity of6.6x 10-"^ ergs"'. Extrap- 
olating this power law model for the nuclear emission to the 
0.5-2 ke V band, we find that the observed count rate exceeds 
that predicted by the model by < 2.4 ct. If we assume this 
excess 0.5-2keV emission is due entirely to a 0.5 - Ike V 
hot gas corona with solar abundance that uniformly fills a 
spherical volume of radius 150 pc, we find an upper limit 
on the mean electron density and gas mass of the gas corona 
of ~ 0.1 cm"^ and < 3 x 10'^ Mq, respectively. These should 
be considered upper limits, because if the nuclear spectrum 
is steeper, the nuclear contribution to the 0.5-2keV emis- 
sion increases. For F > 1.7, the nuclear point source alone 
is sufficient to explain all of the observed soft emission from 
the central region. Thus most of the hot gas associated with 
NGC5846A resides in the stripped tail. 

6. CONCLUSIONS 

NGC 5846 may be the best nearby example of a galaxy 
group in which we can simultaneously study the dynamics of 
non-hydrostatic gas motions induced by galaxy interactions, 
AGN activity and bubble evolution. In this paper we have 
used the combined 120ks Chandra exposure for NGC 5846 
to analyze the signature X-ray surface brightness edges and 
cavities produced by this activity, and have compared these 
data with radio and Ha images of the same region. We also 
observed for the first time ram pressure stripping of a compact 
elliptical galaxy NGC5846A, during its supersonic infall to- 
wards NGC 5846. We find: 

• The surface brightness distribution for NGC 5846 is 
asymmetrically displaced northeast from the nucleus, 
displaying surface brightness discontinuities (edges) at 
19.8 kpc and 6.7 kpc from the nucleus to the northeast 
and at 19.1 kpc and 1 1 kpc to the southwest. Assuming 
a double power law model for the electron density and 
no strong abundance gradients, the density jump at the 
outer (inner) edges are 2.9;!;y 4 (1 .4 ± 0.2) for the edges 
to the northeast and 2.I+0 4 (1.6 ±0.2) to the south- 
east. These multiple surface brightness edges are likely 
the signature of non-hydrostatic gas motions (sloshing) 
caused by the off-axis encounter of NGC 5846 with an- 
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other group galaxy, possibly the large spiral NGC 5850, 
or a merging subgroup. 

• Temperature maps for NGC 5846 confirm that the ob- 
served surface brightness edges correspond to cold 
fronts. Spectral modeling with a VAPEC thermal 
plasma model across the northeast and southwest outer 
edges give temperature ratios of ~ 0.7 and 0.9, respec- 
tively, and pressure ratios ^ 2 suggesting that the gas 
in the cold front may be moving transonically relative 
to its surroundings. Thermal pressures across the inner 
edges are consistent with pressure equiUbrium and zero 
relative velocity. 

• The surface brightness map of NGC 5846 shows a 
spiral-like tail extending 38.5 kpc to the northwest. 
Such spiral structures are also characteristic of gas 
sloshing. However, no corresponding spiral structure 
is observed in the temperature map. This suggests that 
the perturbing encounter did not occur in the plane of 
the sky, and may be consistent with the suggestion by 
Higdon et al. (1998) that NGC 5846 recently experi- 
enced a close off-axis interaction with the spiral galaxy 
NGC 5850. 

• Evidence for recent AGN activity in NGC 5846 are 
seen in two roughly symmetrical bubbles with centers 
~ 0.6 kpc from the nucleus, and in a ghost X-ray cav- 
ity 5.2 kpc to the east. The inner bubbles are coincident 
with 1 .45 GHz and 5 GHz radio emission. Hq emission 
along the inner edge of the northern bubble rim sug- 
gests that gas is cooling, supporting the complex, multi- 
phase nature of gas in the central region. However, 
the Ha emission is weak 6.6 x 10"'"* erg s ' cm"^ ), 
such that the inferred star formation rate is only ~ 
0.037 Moyr-i. 

• From the bubble morphology and thermodynamic prop- 
erties of the surrounding gas, we infer outburst ener- 
gies, bubble age, and AGN mechanical luminosity of 
~ 10^^ ergs, - 2Myr, and (2-4) x lO'^'ergs"' , re- 
spectively, for the inner bubbles, consistent with previ- 
ous results. The outburst energy and age of the ghost 
cavity are ^ 5 x 10^^ ergs and 12Myr, implying an 
AGN duty cycle of ~ lOMyr. Both the inner bubble 
and ghost cavity ages are shorter than the time since 
NGC5850's closest approach to NGC 5846, such that 
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non-hydrostatic bulk gas motions from the encounter 
may have influenced AGN activity, bubble morphol- 
ogy and position, and possibly bubble disruption in 
NGC 5846. 

• The inner bubble rims are threaded with ~ 9 bright, 
knot-like features with individual 0.5-2keV luminosi- 
ties < 6 X lO^** ergs"' . The total 0.5 - 2keV X-ray 
luminosity of these 'knots' is 0.3 x 10'"' ergs"' , com- 
pared to X-ray luminosities of 1.0 x 10"*" erg s"' and 
0.65 X lO'^'ergs"' for diffuse gas in the north and south 
rims, respectively. The mean temperature derived from 
modeUng the stacked spectrum of the 'knots' 

is statistically indistinguishable from that of the sur- 
rounding gas in the rim. The electron density and ther- 
mal pressure in the brightest 'knot', derived assuming 
the 'knot' is spherical, is at least a factor three higher 
than that in the surrounding rim gas. We suggest either 
the 'knots' are extended, filamentary structures in equi- 
librium with the surrounding gas or they are transient 
structures, i.e. gas clouds heated and compressed by 
the recent passage of a shock from the last AGN out- 
burst. 

• Chandra images of the compact elUptical galaxy 
NGC5846A, infalling supersonically towards 
NGC 5846 from the south, show a compact (< 150pc 
radius) corona at the galaxy center and a 0.5 kpc X-ray 
tail containing ~ 10^ of gas, suggesting ongoing 
ram-pressure stripping. Assuming that the observed 
gas originated from passively evolving stars, the time 
between stripping episodes needed to replenish the 
galaxy ISM with sufficient mass to explain the gas 
mass in the tail is ^ 5 .5 Myr. 
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AZIMUTHALLY AVERAGED SPECTRAL FITS 







kT 




null 


(arcsec, arcsec) 


(kpc) 


(keV) 






0,6 


0.35 


72+* ''2 

n 66+" •'2 


76/48 


0.006 


6, 13.2 


1.13 


143/66 


10-^ 


13.2,21.8 


2.03 


n £4+0.02 
U.04 „ „, 

62+^'''^ 
-»■!!? 

h3 

"•°-'-0.02 


72/66 


0.29 


21.8, 32.2 
32.2, 44.6 


3.13 
4.45 


64/67 
56/68 


0.58 
0.86 


44.6, 59.6 


6.04 


65+<'"2 


92/72 


0.06 


59.6, 77.5 


7.95 


70+'' "' 
-Bill 


105/74 


0.01 


77.5, 99.0 


10.24 


73+'' ''' 


89/78 


0.18 


99.0, 124.8 


12.98 


75+'' "i 


81.7/82 


0.49 


124.8, 155.8 


16.27 


n 1+0.02 


110/82 


0.02 


155.8, 192.9 


20.22 


1 01+''''^ 


100/84 


0.11 


192.9, 237.5 


24.96 


I 15+0.06 


113/85 


0.02 


237.5,291.0 


30.65 


1 24+8:M 

^•^^-0.08 


74.6/95 


0.94 


291.0, 355.2 


37.48 


^•^^-0.11 


89/99 


0.75 



Note. — Azimuthally average spectra for the 
NGC 5846 group from Chandra OBSID 7923, taken in 
circular annuli centered on the nucleus of NGC 5846 and 
with logarithmically increasing bin width, rj, ro and Tmean 
are the inner, outer, and mean radius of each annulus. 
All spectra,except the central bin, are modeled with an 
absorbed VAPEC model with fixed Galactic absorption 
(4.24 X 10-" cm"^ ), and abundances for O, Mg, Si.and Fe 
of 0.24, 0.53, 0.72, and 0.492©, re spectiv ely. All other 
abundances are fixed at 0.5 Zq (see gTlJ. 
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TABLE 2 

Angular Sector Density Model Fits 



Sector 


ri 
(kpc) 




(kpc) 


Ji 




"mid 


a2 


NE 
SW 


6.6-6.8 
11.0 


J 4+0.2 

1 6+^-^ 


19.8 
19.1 


2 9+0.5 
2 1+^'^ 

^•^-0.4 


-0 5+0" 

-07^:^ 

'-0.4 


-o.6i^:«? 


-1.2+°' 
-1 2+^-^ 

^•^-0.2 



Note. — Density fit parameters for the double edge model defined in eqns. 1 -|4] 
for the density profiles shown in Fig. [8] Column definitions are: (1) angular sector, (2) 
inner edge location, (3) inner density jump, (4) outer edge location (5) outer density 
jump (6) inner, (7) middle, (8) outer slopes of the density power law. The angular 
sectors are centered at the nucleus of NGC5846 and sul^tend the angle from 123° 
to 180° (NE) and from 268° to 338° (SW), measured counterclockwise from west. 
Uncertainties are 90% confidence levels. 



TABLE 3 

Temperatures Across the NE Edges 



Region 


bi.ai 


bo,ao 


kT 


xVdof 




(arcmin, arcmin) 


(arcmin, arcmin) 


(keV) 




NEl-2 


0.62,0.72 


0.74,0.86 


0.64 ±0.05 


15.8/20 


NEl-1 


0.74,0.86 


0.89,1.03 


0.61 ±0.04 


23.1/22 


NEIO 


0.89,1.03 


1.07,1.24 


0.62 ±0.03 


27.4/23 


NEIl 


1.07,1.24 


1.29,1.48 


0.64 ± 0.03 


39.5/30 


NEO-2 


1.03, 1.19 


1.66, 1.91 


0.63 ±0.02 


48.6/58 


NEO-1 


1.66, 1.91 


2.65, 3.05 


0.71 ±0.02 


69.5/66 


NEOO 


2.65, 3.05 


4.24, 4.88 


1.00 ±0.04 


75.4/57 



Note. — Spectral model fits in concentric elhptical annular regions con- 
centric to the bounding ellipse described in Fig. [8] and restricted to the an- 
gular sector from 123° to 180° measured counter clockwise from west. The 
northeast inner edge is at the boundary between NET- 1 and NEIO. The north- 
east outer edge is at the boundary between NEO-1 and NEOO. The spec- 
tral model is an absorbed VAPEC model with Galactic absorption and abun- 
dances fixed as in Tabled 



TABLE 4 

Temperatures Across the SW Edges 



Region 


ri.ro 

(arcmin, arcmin) 


kT 
(keV) 


x'/dof 


SWI-2 
SWI-1 
SWIO 
SWO-1 
SWOO 
SWOl 


0.95,1.24 
1.24,1.61 
1.61,2.10 
2.11, 2.74 
2.74, 3.57 
3.57,4.64 


76-*^ °^ 

'^-0.05 

94+0.05 
"•^^-0.05 

1 03-^°'^ 
1 17+0.10 
^■^'-0.19 
1 014+<''" 


44.5/34 
11.0/19 
26.1/33 
17.6/36 
25.6/27 
35.6/31 



Note. — Spectral model fits in concentric circular 
annular regions centered on NGC5846's nucleus and 
restricted to lie in the angular sector from 268.2° to 
338.0° measured counter clockwise from west. The 
southwest inner edge is at the boundary between S WI- 1 
and SWIO. The southwest outer edge is at the boundary 
between SWO- 1 and SWOO. The spectral model is an 
absorbed VAPEC model with Galactic absorption and 
abundances fixed as in Table[T] 



Edges and Bubbles in NGC 5846 



TABLE 5 
Cold Front Analyses 



Edge redge "i/no Tt/T,, pi/po Mach 



T q+0.5 71+0-5 9 1+0-5 I n+0.2 S9n+'0 
^•^-0.4 ''-"-OS ^-^-OA ^-^'-02 ^^'-'-100 



(kpc) 


NEouter 


19.8 


SWouter 


19.1 


NF- 
-"^^ -'-'inner 


6.7 


S Winner 


11.0 



(kms-') 



2.lf 0.88:0.21 1.8:1.2 0.9^1 510:|« 
l.4f 0.98!; 1.3!; 0.6^; 250! 
1.6^; 0.79!; 1-3!; 0.6^; 260!^ 



Note. — Cold front ratio analyses following Vikhlinin et al. (2001). 
The edge positions are measured from the center of NGC 5846. Mach 
numbers are relative to the speed of sound outside each edge. 



TABLE 6 

Details of the VLA Radio Observations 



Project 


Observation 


Array 


Frequency 


Bandwidth 


Integration ' 


FWHM, PA 2 


rms 




date 




(GHz) 


(MHz) 


time (min) 


("x", °) 


(MJyb-i) 


AW202 


Jan 1988 


B 


1.5 


50 


12 


5.9 X 4.6, -21 


35 


AF389 


Mar 2002 


A 


1.5 


50 


153.7 


1.5 X 1.5,0 


11.2 


AF142 


Feb 1987 


CnD 


4.9 


50 


3 


8.7 X 4.7, -74 


25 



Note. — ' Integration time on source; ^ full-width half-maximum (FWHM) and position angle (PA) of 
the full array. 



TABLE 7 
Bubble Properties 



Label 


d 


r 


ApV 


ts 


tb 


^mech 




(kpc) 


(kpc) 


( lO'^ ergs) 


(Myr) 


(Myr) 


(lO'*' ergs-l) 


ghost 


5.23 


1.68 


6 


11.2 


12.1 


1.5 


north 


0.75 


0.58 


1 


1.8 


1.1 


3.8 


south 


0.93 


0.58 


1 


2.2 


1.6 


2.5 



Note. — Columns (1) bubble label, (2) distance d from the nucleus of 
NGC 5846, (3) bubble radius, (4) the work needed to evacuate the cavity for 
a relativistic plasma (7 = 4/3), (5) bubble age for bubble rising at the speed 
of sound in the ambient gas (6) bubble age for bubble rising buoyantly at 
its terminal velocity, (7) the instantaneous mechanical power of the outburst 
estimated using the buoyancy timescale 



TABLE 8 

Central Spectral Extraction Regions 



Label 


RA,Dec 


r flout, ^out 






(J2000.0) 


(arcsec) (arcsec,arcsec) 


(arcsec,arcsec) 



Knotl 


15 


06 


29.080, +01 


36 


27.09 


1.09 


Knot2 


15 


06 


29.182, +01 


36 


29.12 


1.20 


Knot3 


15 


06 


29.007. +01 


36 


21.06 


1.06 


Knot4 


15 


06 


28.953. +01 


36 


24.39 


1.02 


Knots 


15 


06 


28.850. +01 


36 


18.90 


1.07 


Knot6 


15 


06 


29.505. +01 


36 


18.86 


1.40 


Knot? 


15 


06 


29.473. +01 


36 


30.65 


1.20 


Knots 


15 


06 


29.329, +01 


36 


29.79 


1.20 


Knot9 


15 


06 


29.078, +01 


36 


19.29 


1.06 


North Rim 


15 


06 


29.518,+01 


36 


24.75 




South Rim 


15 


06 


28.943, +01 


36 


13.93 




Background 


15 


06 


29.223, +01 


36 


21.06 





10.3, 9.6 
8.8, 8.3 
33.2,31.5 



4.5,4.1 
3.7, 3.5 
19.5, 18.4 



Note. — Spectral regions for features in the central region of NGC 5846. Column (1) region 
label, (2) spectral extraction region center, (3) radius of circular regions used for knots and the 
nucleus, (4) outer and (5) inner semi-major, semi-minor axes of elliptical annulus regions used for 
the north and south rims and common local background. Knots 1,2,4,6 — 8 are located only in 
the North Rim, knot 5 is only in the South Rim, and knots 3 and 9 are in the projected overlapping 
region common to both rims. 
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TABLE 9 

Central Features Source Counts 



Label 0.6-0.8keV 0.8-1. IkeV l.l-2keV 2-4keV 





(net counts) 


(net counts) 


(net counts) 


(net counts) 


Knotl 


25±5 


64±9 


25 ±5 


3±2 


Knot2 


24±5 


53 ±8 


44±7 


4±2 


Knot3 


I8±5 


44±7 


28±6 


4±2 


Knot4 


26±6 


43±7 


43 ±7 


4±2 


Knots 


22±5 


36±7 


22±5 


2±1 


Knot6 


39±7 


85 ±10 


52 ±8 


4±2 


Knot? 


21±5 


43±7 


33±6 


-0.37 ±0.03 


Knot8 


18±5 


58±8 


25 ±5 


1±1 


Knot9 


29±6 


56±8 


28 ±6 


14±4 


North Rim 


607 ±29 


1430 ±44 


884 ±33 


113±12 


South Rim 


394 ±24 


922 ±37 


505 ±26 


89 ±10 



Note. — Net source counts in various energy bands for the features 
identified in Table[8]suinmed over both observations, and using the common 
local background also given in Table[8] Knots 1,2, 4, 6 — 8 are located only 
in the North Rim, knot 5 is only in the South Rim, and knots 3 and 9 are in 
the projected overlapping region common to both rims. 



TABLE 10 
Central Features Hardness Ratios 



Label Total Counts (FeS-FeSS)/(FeS+FeSS) (M-S)/(S+M + H) 



Knotl 


117 


0.44 


-0.56 


Knot2 


125 


0.37 


-0.27 


Knot3 


94 


0.43 


-0.36 


Knot4 


116 


0.24 


-0.23 


Knot5 


82 


0.23 


-0.44 


Knot6 


181 


0.37 


-0.40 


Knot7 


97 


0.35 


-0.32 


Knot8 


102 


0.52 


-0.50 


Knot9 


127 


0.32 


-0.44 


North Rim 


3034 


0.40 


-0.38 


South Rim 


1910 


0.40 


-0.42 



Note. — Hardness ratios in the Fe peak and medium energy band for the 
features identified in Tablejs] using the net source counts summed over both ob- 
servations, given in Table |9] FeSS, FeS, M, H and S denote the 0.6-0.8keV, 
0.8-1.1 keV, l.l-2.0keV, 2.0-4.0 keV and 0.6- 1.1 keV energy bands, respec- 
tively, and total counts are the net source counts in the 0.6-4keV energy band. 



TABLE 1 1 
Central Spectral Models 



Label 


kT 
(keV) 


xVdof 


(lO^Ogrgs-i ) 


North Rim 
South Rim 
Stacked Knots 
Knot6 


72+«"3 
64+''"2 

0.72f 


111/65 
26.9/39 
67.4/45 
8.0/15 


1.0 
0.65 
0.31 
0.057 



Note. — Spectral models fit over the 0.5 - 2 keV en- 
ergy range using an absorbed VAPEC mo del wi th Galac- 
tic absorption and abundances fixed as in j|3.2.2l Lx is the 
0.5 — 2keV intrinsic X-ray luminosity, assuming a lumi- 
nosity distance of 24.2 Mpc. Spectra for all regions other 
than the south rim are extracted from OBSID 7923 data 
only. The south rim is a simultaneous fit to OBSID 7923 
and 788. 



